arXiv:1501.00581v2 [astro-ph.GA] 28 Jul 2015 


Draft version July 29, 2015 

Preprint typeset using style emulateapj v. 5/2/11 


THE NATURE AND ORBIT OF THE OPHIUCHUS STREAM 

Branimir Sesar^'^, Jo Bovy®’*, Edouard J. Bernard®, Nelson Caldwell®, Judith G. Cohen^, Morgan Fouesneau®, 
Christian I. Johnson®, Melissa Ness®, Annette M. N. Ferguson®, Nicolas F. Martin*-®, Adrian 
M. Price-Whelan®, Hans-Walter Rix®, Edward F. Schlafly®, William S. Burgett®®, Kenneth C. Chambers®®, 
Heather Flewelling®®, Klaus W. Hodapp®®, Nick Kaiser®®, Eugene A. Magnier®®, Imants Platais®®, John 
L. Tonry®®, Christopher Waters®®, and Rosemary F. G. Wyse®® 

Draft version July 29, 2015 

ABSTRACT 

The Ophiuchus stream is a recently discovered stellar tidal stream in the Milky Way. We present 
high-quality spectroscopic data for 14 stream member stars obtained using the Keck and MMT tele¬ 
scopes. We confirm the stream as a fast moving {vios 290 km s“®), kinematically cold group 
{^vios ^ s“®) of a-enhanced and metal-poor stars ([a/Fe] ^ 0.4 dex, [Fe/H] ^ —2.0 dex). Using a 

probabilistic technique, we model the stream simultaneously in line-of-sight velocity, color-magnitude, 
coordinate, and proper motion space, and so determine its distribution in 6D phase-space. We find 
that that the stream extends in distance from 7.5 to 9 kpc from the Sun; it is 50 times longer than 
wide, merely appearing highly foreshortened in projection. The analysis of the stellar population 
contained in the stream suggests that it is ~ 12 Gyr old, and that its initial stellar mass was ~ 2 x 10"® 

Mq (or at least > 7 x 10^ Mq). Assuming a fiducial Milky Way potential, we fit an orbit to the 
stream which matches the observed phase-space distribution, except for some tension in the proper 
motions: the stream has an orbital period of ^ 350 Myr, and is on a fairly eccentric orbit (e ~ 0.66) 
with a pericenter of ^ 3.5 kpc and an apocenter of ^ 17 kpc. The phase-space structure and stellar 
population of the stream show that its progenitor must have been a globular cluster that was dis¬ 
rupted only ^ 240 Myr ago. We do not detect any significant overdensity of stars along the stream 
that would indicate the presence of a progenitor, and conclude that the stream is all that is left of the 
progenitor. 

Subject headings: globular clusters: general — Galaxy: halo — Galaxy: kinematics and dynamics — 

Galaxy: structure 


1. INTRODUCTION 

One of the main goals of Galactic astronomy is the 
measurement of the Milky Way’s gravitational potential, 
because knowledge of it is required in any study of the 
dynamics or evolution of the Galaxy. An important tool 
in this undertaking are stellar tidal streams, remnants 
of accreted Milky Way satellites that were disrupted by 
tidal forces and stretched into filaments as they orbited 
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in the Galaxy’s potential. The orbit of a stream is sen¬ 
sitive to the properties of the potential and thus can be 


spanned by the stream (e.g 

Koposov et al. 112010 

New- 

berg et al.j2010 Sesar et al. 

2U13| Belokurov et al. 

2014|). 


In this c ontext, the recently d iscovered Ophiuchus stellar 

stream (Bernard et al.|2014b) is very interesting because 


it is located fairly close to the Galactic center (galacto- 
centric distance of ~ 5 kpc), and as such probes the part 
of the potential that other known stellar tidal streams do 
not probe. 

The Ophiuchus stream is a ~ 2.5° long and 7' wide 


stellar stream that was recently discovered by Bernard 


et al. ( 2014b) in the Pan- STARRSl photometric catalog 
(FSl; Kaiser et al. 2010). Bernard et al. inferred from 
its color-magnitude diagram (GMD) that it is consistent 
with an old (> 10 Gyr) and relatively metal-poor pop¬ 
ulation {[Fe/H] ^ —1.3 dex) located ^ 9 kpc away at 
(/, b) ~ (5°, -(-32°). They did not detect a progenitor (or 
a remnant of it), but suggested that the progenitor would 
most likely be a globular cluster. 

Due to the lack of proper motion and line ofsight ve¬ 
locity measurements, Bernard et al. could not determine 
the orbit of stream and thus could not use it to constrain 
the potential. Furthermore, without knowing the orbit 
of the stream, they could not fully explain two interest¬ 
ing properties of the Ophiuchus stream, namely, its very 
short length and the lack of a visible progenitor. The 
projected angular length of 2.5° at a distance of ^ 9 kpc 
implies a projected physical length of ^ 400 pc for the 
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Ophiuchus stream. Such a short length suggests that 
its progenitor must have been disrupted fairly recently. 
However, if that was the case, the progenitor should still 
be visible as an overdensity of stars somewhere along the 
stream. Yet, no progenitor has been detected so far. 

To address the above questions, we need to know the 
orbit of the Ophiuchus stream, and to determine its or¬ 
bit we need to measure the stream’s line of sight ve¬ 
locity, distance, and proper motion. In Section we 
describe the data we use in this work; the PSl photom¬ 
etry and astrometry, the spectroscopic follow-up of can¬ 
didate stream members, and the measurement of their 
line of sight velocities, element abundances, and proper 
motions. In Section]^ we provide a detailed characteri¬ 
zation of the stream in position, velocity, and abundance 
(7D) phase space. The constraints obtained in Section]^ 
are then used to constrain and examine the orbit of the 
stream (Section® and the time of disruption of its pro¬ 
genitor (Section®. In Section® we discuss the nature of 
the stream’s pecmiar orbit, hi^light the solved and un¬ 
covered puzzles related to the stream, and finally present 
our conclusions. 


2. DATA 


2.1. Overview of the Pan-STARRSl survey 

The PS I survey has observed the entire sky north of 
declination —30° in five filters coverin g 400 — 1000 nm 


Stubbs et al. 

2010 

Tonry et al. 

L '2 - r 

2012 


billion-pixel camera (IHodapp et al. 

2004| lOnaka et al. 

2008| Tonry & Onaka|20(15 ). In singl 

e-epocn images, the 


ratio (S/N) of 5 at 22.0, 22.0, 21.9, 21.0, and 19.8 mag 
in PSl grizypi bands, respectively. The survey pipeline 
automatically processes images and perf orms photome- 


try and astrometry on detected sources ( 

Magnier 

2006 

20071. The uncertainty 

in photometric cat 

ibration 

ot the 

le as- 

survey is ^ 0.01 mag ( 

Schlafly et al. 

2012 

), and t 


trometri c precision of ^gle- epoch detections is 10 mil- 


liarcsec (Magnier et al.|[2008 hereafter mas). 


2.2. Line of sight velocities 


Based on the findings of Bernard et al. (2014b), we 
have used the deredd ened fiducial of the old globular clus¬ 
ter NGC 5904 (from Bernard et al.|2014a and shifted to 
the distance modulus of 14.9 mag) to select ~ 170 can¬ 
didate Ophiuchus stream members from the PSl pho¬ 
tometric catalog. The candidates were selected if their 
position was within 4.5' of the best-fitt ing great circle 
contain ing the stream (see Figure 1 of Bernard et al. 
2014b), and if their dereddened gpi — ipi color and ipi- 
band magnitude were within 0.1 mag and 0.5 mag of 
the fiducial isochrone, respectively. We observed the se¬ 
lected ca ndidates using the DEIMOS spectrograph on 
Keck B (Faber et al. 2003) and using the Hectochelle 
fiber spectrograph on MMT (Szentgyorgyi et al. 2011) 
over a course of two nights. 

Seven candidate blue horizontal branch stars were ob¬ 
served with DEIMOS on 2014 May 29th (project ID 
2014A-C171D, PI: J. Cohen). The observations were 
made using the 0.8" slit and the high resolution (1200 
G) grating, delivering a resolution of 1.2 A in the 6250- 
8900 A range. The spectra were extracted and calibrated 


using standard IRAlj^ tasks. The uncertainty in the 
zero-point of wavelength calibration (measured using sky 
lines) was < 0.04 A (< 2 km s“^ at 6563 A). 

The line of sight velocities of stars observed by 
DEIMOS were measured by fitting observed spec tra with 
synthetic te mplate spectra selected from the |Munari| 


et al. (2005) spectral librarjJ^ Prior to fitting, the syn- 
thetic spectra were resampled to the same A per pixel 
scale as the observed spectrum and convolved with an 
appropriate Line Spread Function. The velocity obtained 
from the best-fit template was corrected to the barycen- 
tric system and adopted as the line of sight velocity, uios- 
We added in quadrature the uncertainty in the zero-point 
of wavelength calibration (2 km s“^ at 6563 A) to the ve¬ 
locity error from fitting. 

The remaining 163 Ophiuchus stream candidates were 
observed with Hectochelle on 2014 June 6th (proposal ID 
2014B-SAO-4, PI: G. Johnson). Observations were made 
using the RV31 radial velocity filter, which includes Mg 
I/Mgb features in the 5150-5300 A range. To improve 
the S/N of faint targets, we binned the detector by 3 
pixels in the spectral direction, resulting in an effective 
resolution of i? ~ 38,000. 

Hect ochelle spectra were extracted and calibrated fol¬ 
lowing Caldwell et al. (2009). To account for variations in 
the fiber throughput, the spectra were normalized before 
sky subtraction. The normalization factor was estimated 
using the strength of several night sky emission lines in 
the appropriate order. Sky subtraction was performed 
using the average of 20-30 sky fi bers, using the method 
devised by Koposov et al. (2011). A comparison of ob¬ 
served and laboratory positions of sky emission lines did 
not reveal any significant offsets in wavelength calibra¬ 
tion (i.e., no offsets greater than 0.5 km s“^ at 5225 A). 

The line of sight velocities of stars observed b y Hec¬ 
tochelle were measured using the RVSAO package (|Kurtz| 
'fc Mink|1998 ), by cross-correlating observed spectra with 


a synthetic spectru m of an A-type and a G-type giant 


star (constructed by Latham et al.||2002|). The velocity 
obtained from the best fitting template was adopted. To 
the uncertainty in uios, measured by RVSAO, we added 
(in quadrature) the uncertainty in the zero-point of wave¬ 
length calibration, which we measured using sky emission 
lines to be 0.5 km s“^. Finally, the measured velocities 
were corrected to the barycentric system using the BCV- 
CORR task. 

A comparison of velocities measured from DEIMOS 
and Hectochelle spectra for star “bhb6” (Table |^, shows 
that the two velocity sets are consistent within stated 
uncertainties. 


2.3. Element abundances 

Even though the primary goal of spectroscopic obser¬ 
vations was to obtain precise radial velocities, the wave¬ 
length range and the resolution of Hectochelle spectra 
are sufficient to allow estimates of element abundances. 

We determined stellar parameters from the continuum- 
normalized, radial velocit y-corrected spectra using the 
' Gasevl(l2014D, 
ienl (119731) . Kr 


SMH code of 


code of Sneden 11973)' 


|, which is built on the MOOG 
Curucz a—enhanced (0.4 dex) 


http://iraf.noao.edu/ 

http://archives.pd.astro.it/2500-10500/ 
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model atmospheres (Castelli fc Kurucz 


2004) and a 


line li st compiled from h'rebei et ai. ” 120111)07 Casey 
(2014) were used (see Table 5 in the electronic version 
ot the Journal). First, effective temperatures were cal¬ 
culated from 2 MASS photometry and color-temperature 


calibrations of Gonzalez Hernandez &: Bonifacio (2009) 


and spectroscopic temperatures were optimized around 
this value using the SMH code, by removing abundance 
trends with line excitation potential. Parameters of log g 
and [Fe/H] were determined using 14 Fe I and 2-3 Fe II 
lines to achieve ionization balance, and microturbulence 
was calculated by removing trends in abundances as a 
function of the reduced equivalent width of the lines. Es¬ 
timates of the a—enhancement were obtained using only 
the few available clean Mg, Ca and Ti (I and II) lines, 
which comprised a total of 6-8 lines per star. The abun¬ 
dances were measured from equivalent widths and the 
lines we used are listed in Table 5. Because Mg lines are 
strong and may be saturated, the values of [Mg/Fe] are 
significantly different than values of [Ca/Fe] and [Ti/Fe]. 



Figure 1. This plot illustrates the systematic offset in declina¬ 
tion (ADecl) of objects observed in the POSS-II Blue epoch of 
the USNO-B catalog (plate 799), as a function of the ^pi-band 
magnitude. For clarity, only a subset of objects are plotted. The 
solid yellow circles show the median AD eel in magnitude bins, and 
the solid red circles show the rms scatter in magnitude bins. Note 
how the brighter objects are systematically offset by ~ 200 mas 
from the fainter objects. The rms scatter indicates that the aver¬ 
age precision in this coordinate and epoch is ~ 120 mas for objects 
brighter than gpi = 19. 


2.4. Proper motions 

Proper motions are crucial constraints for determ ining 
the orbit of a stream (e.g., |Koposov et al.||2010|). To 
measure the proper motion of stars in the vicinity of the 
Ophiuchus st ream, we combine the astrometr y provided 


by USNO-B (|Monet et al. 2003) and 2MASS (Skrutskie 
et al.|2006 l catalogs with the PSl catalog. The USJNU-B 
catalog lists photometry and astrometry measured from 
photographic plates in five different band-passes (O, E, 
J, F, and N). The USNO-B plates were exposed at 
different epochs, and thus each object in the catalog can 
have a maximum of five recorded positions. The 2MASS 
catalog provides only one position entry per object. 

To reduce the systematic offsets in astrometry between 
different catalogs, we first calibrate USNO-B and 2MASS 
positions to a reference frame that is defined by positions 
of galaxies observed in PSl (which are on the ICRS coor¬ 
dinate system). We define galaxies as objects that have 
the difference between point-spread function (PSF) and 
aperture magnitudes in PSl rpi and ipi bands between 
0.3 and 1.0 mag. 

The astrometric reference catalog is created by aver¬ 
aging out repeatedly observed positions of PSl objects. 
Between 2012 May and June, the region in the vicinity 
of the Ophiuchus stream was observed four times in PSl 
ffpi, ?'pi) and ipi bands. To minimize the uncertainty 
in astrometry due to wavelength-dependent effects, such 
as the differential chromatic refraction, we only average 
out positions observed through the rpi-band filter. Since 
the a strometric precision o f single-epoch detections is 10 
mas (Magnier et al. 2008), the precision of the average 

position is ^ 5 mas or better. _ 

The USNO -B astrometry is calibrated following |Munn| 


et al. 
late t 


(2004 see their Section 2.1). First, we calcu- 
le positions of objects at each of the five USNO- 


B epochs, using software kindly provided by J. Munn. 
Then, for each USNO-B object we find the 100 nearest 
galaxies, calculate the median offsets in right ascension 
and declination between the reference PSl position and 
the USNO-B position for these galaxies, and add the off¬ 
sets to the USNO-B position in question. This is done 
separately for each of the five USNO-B epochs. The 
single-epoch 2MASS positions are calibrated using the 


same procedure. 

Having tied the positions for each object at one 2MASS 
and five USNO-B epochs to the PSl astrometric reference 
frame, we can now check for any additional systematic 
uncertainties in the calibrated astrometry. We do so us¬ 
ing the leave-one-out cross-validation. One of the six 
calibrated positions is withheld, and a straight line is fit¬ 
ted to the remaining five positions and the PSl position. 
The straight line fit (i.e., essentially a proper motion fit, 
neglecting the parallax) is then used to predict the posi¬ 
tion of an object at the withheld epoch. The difference 
between the withheld position and the predicted position 
is labeled as Ai?A or ADecl. 

Inspection of ARA or ADecl values has revealed that 
the positions of USNO-B objects depend on magnitude 
for some epochs (Figure]^. We have examined ARA and 
ADecl values in different regions of the sky, and have 
concluded that these astrometric issues affect individual 
photographic plates, and are not specific to a particular 
photographic bandpass. To remove this dependence, we 
subtract plate-specific and magnitude-dependent offsets 
(shown by yellow circles in Figure from the original 
USNO-B positions before we calibrate the positions us¬ 
ing PSl galaxies. Since USNO-B does not provide un¬ 
certainty in positions, we adopt the rms scatter of ARA 
and ADecl values (shown by red solid circles in Figure[^ 
as an estimate of the uncertainty in position at a given 
magnitude and epoch. 

We have also examined whether ARA and ADecl as¬ 
trometric residuals depend on the gpi — ipi color. We 
find that the residuals do depend on the color, and that 
they can be as high as 100 mas. The most likely expla¬ 
nation for this dependence is the differential chromatic 
refracation. We correct for this dependence using a sim¬ 
ilar approach as above. For each photographic plate we 
bin ARA and ADecl residuals as a function of color, cal¬ 
culate the median for each color bin, and subtract that 
value from the observed positions of stars in that color 
bin, before we calibrate the positions using PSl galaxies. 

Finally, to measure the absolute proper motion of an 
object we fit a straight line to all available positions 
(min. 3, max. 7) as a function of time (the epoch baseline 
is 58 years). The proper motion of confirmed Ophiuchus 
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Figure 2. The distribution of heliocentric line of sight velocities 
of stars observed by DEIMOS and Hectochelle. The uncertainty 
in individual vio^ measurements is < 2 km s~^ and the bin size is 
25 km s“^. The Ophiuchus stream is detected as a group of stars 
with vios ~ 290 km s“^. 


stream members is listed in Table U) 

For verification, we have also measured the proper 
motion of 700 cand idate QSOs, selected using WISE 
( Wright et al. W1 — W2 > 0 .8 color criterion (see 

Section 2.1 by JNikutta et al.||2014 ). The median proper 
motion of candidate QSOs is 0.3 mas yr“^ and the uncer¬ 
tainty of the median is 0.2 mas yr“^, showing that there 
is no statistically significant offset in measured proper 
motions. 

To measure the systematic uncertainty in proper mo¬ 
tions, one should ideally calculate the rms scatter of 
proper motions of fairly bright and static sources. Bright 
QSOs would be an ideal choice for this measurement, be¬ 
cause they are extragalactic point sources and because 
they are not used in the calibration process. Unfortu¬ 
nately, the candidate QSOs described above are too faint 
(r > 18 mag) to be used for this purpose (i.e., the uncer¬ 
tainty in their proper motions is already dominated by 
Poisson noise). 

Instead, we measure the systematic uncertainty by 
calculating the median uncertainty in proper motion of 
bright stars. For stars brighter than r = 17 mag, the un¬ 
certainty in proper motion is ~ 1.5 mas/yr (and constant 
with magnitude), and we adopt this value as the system¬ 
atic uncertainty. To determine whether this uncertainty 
is well-measured, we examined the distribution of per- 
degrees of freedom values calculated from proper motion 
fits of bright stars (r < 17). The mode of this distri¬ 
bution is centered at 1, indicating that on average, the 
uncertainties in proper motion are well-measured (i.e., 
not overestimated or underestimated). However, the ob¬ 
served distribution has longer tails than expected (to¬ 
ward high values), indicating that for some stars the 
uncertainties in proper motion are understimated (e.g., 
due to blending of sources in photographic plates). 


3. CHARACTERIZATION OF THE OPHIUCHUS STREAM 
3.1. Line of sight velocities 

The vios distribution of stars observed by DEIMOS and 
Hectochelle is shown in Figurej^ In this figure, a group of 
14 stars with 285 < z;;os/kins“^< 292 clearly stands out. 
This group, which we identify as the Ophiuchus stream, 
is well-separated from the majority of stars which have 
|r!ios| < 200 km s“^. The positions, velocities, and PSl 
photometry of stars in this group are listed in Table 
A closer look at vios of stars in the Ophiuchus stream 


Figure 3. Line-of-sight velocities of stars in the Ophiuchus 
stream are shown as symbols with error bars. The thick solid line 
shows the most probable model {vios{^) = 4.0(£ — 5) + 289.1 km 
s~^). To illustrate the uncertainty in the most probable model, the 
thin semi-transparent red lines show 200 models sampled from the 
posterior distribution. The vert ical dashed lines show the likely 
extent of the stream (see Section [3^. 


(Figure]^ suggests that their velocities are changing as a 
function of galactic longitude. To fit this possible velocity 
gradient, we use an a pproach similar to the one taken by 


Martin & Jin (2010 see their Section 2.1.1). 

We wish to nnd the most likely set of parameters 0 for 
which the observations of stars listed in Table I? = 
{dk}i<fc<i 4 , match the model described below. In the 
current problem, each data point dk is defined by its 
line of sight velocity vios,k at galactic longitude (.k^ dk = 
{vios,k,^k}- The velocity has an uncertainty of 
The star “bhb6” has been observed twice, so for this 
star we adopt the weighted average (and its associated 
uncertainty) of the two line of sight velocity observations. 
The uncertainty in longitude is not considered because it 
is much smaller compared to the uncertainty in velocity. 
The data points are also considered to be independent. 
Therefore, the likelihood of these data points with the 
model defined by the set of parameters 9, is 


p{v\e) = l[p{d^\e,£^), (1) 

k 


where p(dk|0,I'k) is the likelihood of data point k to be 
generated from the model. Using Bayes theorem, the 
probability of a model given the data, p{9\'D), is 

pie\V)^p{V\0)p{9), (2) 

where p{9) represents our prior knowledge on the model. 
We explicitly define the likelihood p(dk|0,€k) as 

p{i^los,k\9^ — Ifi^Vios ^k') ^ (^) 

where 

M{x\p, a) = (l/vT^) exp(—0.5((a: — ^)^)/cr^) (4) 

is a normal distribution, and 9 = are pa¬ 

rameters of our model: 


1. is the velocity gradient along the galactic lon¬ 
gitude direction, 

2. Vios is the velocity of the stream at the reference 
galactic longitude £o = 5°, and 

3. s accounts for the additional scatter in velocities 
(e.g., due to the intrinsic velocity dispersion of the 
stream). 
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The f(ffe) is the predicted velocity of the stream at po¬ 
sition £k 

v{ik) = -^(4 - 4) -\-W^ (5) 

at 


-I- cr^ is the quadratic sum of the addi- 


and cr^ = y' 

tional scatter s in velocity and the uncertainty in line of 
sight velocity of data point k. The likelihood of all data 
points can be calculated using Equation 
Before we can sample from the posterior probabil¬ 
ity distribution over our model parameters, we need 
to define the prior probabilities of model parameters. 
As prior probabilities, we adopt priors that are uni¬ 
form in these ranges: 270 < ul77/kms“^ < 320, —8 < 
%^/kms-Meg-i < 8, 0 < s/kms-i < 3. 

T o efficiently explore the p arameter space, we use 
the Goodman & Weare (2010| Affine Invariant Markov 
chain Monte Carlo (iVlCivlC) Ensemble sampler as imple- 
mented in the emcee packagq^ (v2.1, Eoreman-Mackey 


et al. 20131. We use 200 walkers and obtain conver- 


genc^^^l aft^er a short burn-in phase of 100 steps per 
walker. The chains are then restarted around the best-fit 
value and evolved for another 2000 steps. To enable easy 


reconstruction of the posterior distribution, we provide 

chains in a Zenodo data repositorjp® ( 

Sesar et al. 

2015 

doi: 10.528Uzenodo. 19197) for ail of tl 

tie data modeling 


We characterize the most probable set of parameters 
for our model using the maximum a posteriori values. 
We also report the median and equivalent l-cr confidence 
intervals using the 50th, 16th and 84th percentiles, re¬ 
spectively (see Table [^. 

We find = 4.0±1.2 km s“^ deg“^, = 289.1 ± 
0.4 km s“^, and a very small velocity dispersion of s = 
0.4t°;| kms-b Thus, we detect a gradient in the line of 
sight velocity at a 4cr level. 

3.2. Element abundances 

The preliminary element abundances of five red giant 
branch (RGB) stars in the Ophiuchus stream (that is, the 
vios ~ 290 km s“^ group), and observed by Hectochelle, 
are listed in Tabled The uncertainties of the determined 
parameters are listed in the notes of Table 

We Hnd the stars in the Ophiuchus stream to be poor 
in Fe {[Fe/H] ^ —2.0 dex) and enhanced in a-elements 
([of/Fe] = 0.4 ± 0.1 dex). Their [Fe/H] are consistent 
within 0.05 dex (rms scatter), despite fairly large esti¬ 
mated uncertainties in individual measurements (;< 0.2 
dex). The small scatter in [Fe/H] suggests that these 
stars come from the same single stellar population. 

Based on their p osition (within 5' of the Ophiuchus 
stream, as traced by Bernard et al. 2014b), kinematic and 
chemical properties, we conclude that all stars listed in 
Table are high-probability members of the Ophiuchus 
stream. 


3.3. Color-magnitude diagram 


Table 1 

Ophiuchus stream parameters 


Parameter 

MAP“ 

Median and central 68% C.I.^ 

'^los 

289.1 

289.llg'4 km s-l 

d'^los 

d£ 

4.0 

4.0^)'2 km deg“^ 

S 

0.0 

0.4^°;® km s-l 

[Fe/H] 


-1.95lpg dex 

[a/Fe] 


0.4l°'( dex 

Age 

11.7 

Gyr 

Mass-loss parameter rj 

0.49 


Metallicity Z 

2.3 X lO-'* 

(2.3l0;4) X 10-1 


0.011 

0-008l°;°°® mag 

DM 

14.57 

14.58l° °| mag 

dDM 

di 

-0.20 

-0.201oq| mag deg-l 

^g,iso 

0.0003 

0-012l°:°J® mag 

^r,iso 

0.0003 

0-009l“ °o6 mag 

^i,iso 

0.0006 

0-007l°;°°| mag 

^z,iso 

0.0003 

0.006l“;0°| mag 

^y,iso 

0.0003 

0-007l°;°“ mag 

P 

^min 

3.81 

3.84l°;03 deg 

^max 

5.85 

5.86lg;“ deg 

A 

31.37 

31.38tg g2 dgg 

B 

-0.80 

-0.80«°3 

C 

-0.15 

-O.lOtg'Dl deg-i 

Deprojected length 

1.6 

l.5l°;3 kpc 


6.0 

6-9^0 g arcmin 

Ft 

-5.5 

—5.6^q' 3 mas yr“^ 

dm 

~W 

-1.5 

— g mas yr-l deg-l 

Fb 

2.4 

^•^-0 3 

dm 

df 

2.0 

2.3l°;4 mas yr-l deg-l 


Pericenter 

3.55 

3.5Vrohf°o:?^)Apc 

Apocenter 

17.0 

kpc 

Eccentricity 

0.66 

U.OO_Q 01 \^-0.08 J 

Orbital period 

351 

346l).l (tlPj Myr 

Radial period 

239 

237t^ Myr 

Vertical period 

346 

3421)1 (1)5) Myr 

Mass of the progenitor 


~ 2 X lOl Mq 


^ Maximum a posterior value, where available. 

^ The median and the central 68% confidence intervals are measured from 
marginal posterior distributions. The intervals are calculated from the dif¬ 
ference of the 16th and 50th, and 84th and 50th percentile. 

^ Recall that E{B — V)^f f is the offset with respect to the reddening 
provided by the |Schlegel et al.| | [T998^ dust map, and is not reddening by 
itself. 

^ The numbers in parenthesis illustrate the range of values (with respect 
to the median) obtained when varying the distance of the Sun from the 
Galactic center and the circular velocity at solar radius (see Section 

The sample of Ophiuchus stream members, which we 
have identified above using velocities and metallicities, 
now gives us an opportunity to further constrain the dis¬ 
tance and the CMD of the stream. 


http://dan.iel.fm/emcee/current/ 

We checked for convergence of chains by examining the auto¬ 
correlation time of the chains per dimension, 
https://zenodo.org/record/19197 


3.3.1. Model 

To model the CMD of the stream, we use a proba- 
bilis tic a pproach analogous to the one described in Sec¬ 
tion 3.1 In our data set, P, each data point dk is now 
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Table 2 

Element abundances of Ophiuchus stream stars 


Name 

Eff 

(K) 

logs 

(dex) 

[Fe/H] 

(dex) 

[Mg/Fe] 

(dex) 

[Ca/Fe] 

(dex) 

[Ti/Fe] 

(dex) 

rgbl 

5680 

3.2 

-2.00 

-0.12 

0.63 

0.69 

rgb2 

5430 

3.3 

-1.95 

0.06 

0.62 

0.47 

rgb3 

5700 

3.1 

-1.95 

0.23 

0.27 

0.46 

rgb4 

5400 

2.8 

-1.90 

0.24 

0.57 

0.42 

rgb5 

5720 

2.9 

-1.90 

-0.14 

0.83 

- 


Note. — The uncertainty in Te// is < 200 K, < 0.45 dex for log g, 
< 0.2 dex for [Fe/H], and ~ 0.35 dex for abundances of q— elements. 


defined by its galactic longitude and latitude, and by its 
PSl grizypi magnitudes, dk = {ik,bk,gk,rk,ik,Zk,yk}- 

Our data set contains only the Ophiuchus stream stars 
that were identified based on spectroscopic data (i.e., ve¬ 
locity and metallicity). Thus, the set is uncontaminated 
but very sparse and has a complicated spatial selection 
function. Because of this, and because we are primar¬ 
ily interested in constraining the distance of the stream, 
we focus on finding the isochrone(s) that match the con¬ 
firmed members in the color-magnitude space, and do 
not to model the projected s hap e of the stream on the 
sky (for now, but see Section 3.4). 

To model the stream in color-magnitude space, we use 
a grid of theoretical PARSEC isochron eFd (release vl.2S; 
Bressan et al.|[2M2 Chen et al.|2014 ). Each isochrone I 
provides PSl magnitudes m' = g'p^,rp^, i'p^, z'p^.y'p^ for 
a star of initial mass Minu in a single stellar population 
of age t, metal content and parametrized for the mass- 
loss on the RGB using the variable val ues of the Reimers 

119771). At the reference 


law parameter rj ( Reimers| 


galactic longitude £o = S'", the distance modulus of the 
stellar population is defined with parameter DM, and a 
gradient in distance modulus with galactic longitude is 
modeled with parameter . We model the extinction 
in a PSl band by adding the 

Ce^t ■ [Esfd{B - V\tkM) + E{B - V)off] ( 6 ) 
term to isochrone magnitudes, where Espn(B — V\£k, bk) 


IS 


et al. 
tor 


{£k,bk) in the 
E{B - VUf 
offset. The extinction coefli- 


Schlegel 

accounts 


the redde ning at position 
(|1998|) dust map, and 
possible zero-point 

Cea:t = 13.172,2.271,1.682,1.322,1.08 7} for PSl 


a 

cients 

ig,r,i,z,y} _ 

FinkbeinerJ (|2011|). To account for the fact that stel- 


bands are taken from Table 6 of Schlafly & 


lar evolution models are not perfect, we introduce five 
Cm,iso parameters, where m = g,r, i, z, y, that model the 
uncertainty in each PSl grizy magnitude provided by 
PARSEC isochrones. 

Given the above model of the stream, the likelihood of 
a star k with this model is defined by the comparison of 
the spectral energy distribution {g, r, j, z, j/jpi with the 
prediction {g', r',i',z',y'}p^ of an isochrone Ik for a star 
with the initial mass Therefore, 


p(dk|Ife) = J 


A/'(TOfc|TO'(M,™t),cr(„J 


m—g,r,i,z,y 


(7) 


where 

Xk =X{ik,bk,t,Z,r],E{B - V)off,DM, ,aiso) 

is the isochrone at the galactic position of star fc, aiso = 
{o-g,iso,crpiso,a-i,iso,cr^^iso,o-y^iso}, A/'(x|/r, cr) is a normal 
distribution, and 




Cmfc T ^m,iso 


[0.lCextEsFD{B — V\£k,bk)] 

. 

is the sum of the uncertainty in the isochrone magnitude 
(cm.iso), observed magnitude of data point k (cTm;,), and 
extinction (10% frac tional uncertainty in Esfd{B — V)\ 


Schlegel et al. 
can then be ca 


1998). The likelihood of all dat^pointe 
culated by combining Equations and 


3.3.2. Priors on the CMD model 

Before we can calculate the probability of a model, we 
need to define the prior probabilities of model param¬ 
eters. Below, we list our priors and describe the justi¬ 
fication for each one. A summary of priors is given in 
Table I 

Based on spectroscopic data, the Ophiuchus 
stream is metal-poor ([Fe/H] =—1.95 ± 0.05 dex) 
and g—enhanced ([ a/Fe] = 0.4 ± 0.1 dex). As shown by 


Salaris et al. 
models 


(1993 


the a—enhanced stellar population 
are equivalent to scaled-solar ones with the 
global metal content [M/H], where [M/H] for 
models can be calculated using their 


same 

a—enhanced 
Equation 3 

[M/H] « [Ee/i7] -blogio(0.638 x lO'"/-^^! -f 0.362). (10) 

For the element abundance of the Ophiuchus stream, we 
estimate [M/H] ~ —1.7 ± 0.2 dex. This means that we 
should adopt the normal distribution N{\ogiQ{Z/Z q) \ — 
1.7,0.2) as the prior probability of metallicity Z (where 
Zq = 0 .0152 is the solar metal content used by this per- 
ticular set of PARSEC isochrones). However, in the con¬ 
text of cross-validating our analysis, we decided to re¬ 
place the above metallicity prior in favor of a (less infor¬ 
mative) prior that is uniform in the 0.0001 < Z < 0.0004 
range. Even though a less i nformative prior was adopted, 
at the end of Section [3.3.3| we find a very impressive con¬ 
sistency between the posterior distribution of metallicity 
Z (obtained using CMD fitting) and the spectroscopic 
estimate of Z (see bottom panel of Figure [^. In the 
end, it is important to note that our results would not 
have changed significantly if we used the more informa¬ 
tive prior for metallicity content Z. 

The presence of BHB stars, the [Fe/H] and the 
a—enhancement of the stream point to an old stellar pop¬ 
ulation. Thus, for age we adopt a uniform prior in the 
8 < t/Gyr < 13.5 range. 


& Fusi Pecci 

00 

00 

0 ^ 

1 — 1 

eter rj we ac 

opt a 


0.4 (Renzini 


< 77 < 0.5 

range. For the uncertainty in isochrone magnitudes, we 
adopt a prior that is uniform in the 0 < < 0.1 

mag range, whe re m = q,r, i, z, y. 


According to Bernard et al. 


(2014b), the Ophiuchus 


http://stev.oapd.inaf.it/cnid 


stream is located about 9 ± 1 kpc troin the Sun. Thus 
for DM we adopted a uniform prior in the 14.2 < DM < 
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Table 3 

Prior probabilities of CMD parameters 


Parameter 

Prior type 

Range 

Age t 

uniform 

8 to 13.5 Gyr 

Mass-loss parameter r} 

uniform 

0.2 to 0.5 

Metallicity Z 

uniform 

0.0001 to 0.0004 

E{B - VU f 

uniform 

-0.1 to 0.1 mag 

DM 

uniform 

14.2 to 15.2 mag 

dDM 

d£ 

uniform 

-0.5 to 0.5 mag deg~^ 

o-m,iso {m = g,r,i,z,y) 

uniform 

0 to 0.1 mag 


15.2 mag range (corresponding to the 7-11 kpc range). 
For the gradient in distance modulus, a uniform prior 
in the < 0.5 mag deg“^ range is adopted. Since 

the reddening in the region of interest is greater than 0.1 
mag, we assume that the possibl e systematic offset in 
E{B — V) values provided by the iSchlegel et al.| (|1998[ 
dust map is less than 0.1 mag (i.e., \E[B — V)off \ < 0. 
mag). 


3.3.3. Posterior distributions of CMD parameters 

To efhciently explore the parameter space, we again 
use the emcee package. We use 1000 walkers and ob¬ 
tain convergence after a short burn-in phase of 100 steps 
per walker. The chains are then restarted around the 
best-fit value and evolved for another 4000 steps. The 
maximum a posterior values, the median and the central 
68% confidence intervals of model parameters are listed 
in Table [H 

We find the stream to be ^ 12 Gyr old and to have 
a distance modulus of 14.58 ± 0.05 mag (i.e., a distance 
of 8.2 kpc) at the reference galactic longitude Iq = 5° 
(top panel of Figure]^. Most importantly, we detect a 
gradient of —0.20±0.03 mag deg“^ in distance modulus. 
This gradient is inconsistent with zero (i.e., with the no 
gradient h ypothesis) at a 7 a level , and confirms the sug¬ 
gestion by Bernard et al. (2014b|) that the_eastern part 
of the stream is closer to the Sun (Figure]^. The gra¬ 
dient in distance modulus is mostly constrained by BHB 
and MSTO/SGB stars. When these stars are not used to 
constrain the CMD of the stream, the marginal posterior 
distribution of the gradient in distance modulus becomes 
multimodal and poorly constrained. 

To verify whether the observed gradient in distance 
modulus is real, we compared derredened colors and mag¬ 
nitudes of BHB stars in the Ophiuchus stream. Stars 
bhbl and bhb3 have identical dereddened gpi — ipi color 
{gpi — ipi = —0.43 ma g), and thus should have iden¬ 
tical absolute magnitude^^ However, their dereddened 
ipi-band magnitudes differoy 0.1 mag, and the brighter 
star in the pair is located 0.5 deg eastward (in the galac¬ 
tic longitude direction), in agreement with the observed 
distance modulus gradient of -0.2 mag deg^. The stars 
bhb6 and bhb7 also show similar behavior. 

In Figure we compare CMDs of the Ophiuchus 
stream and fiO^d stars. The CMD of field stars (grayscale 
pixels) was obtained by binning the gpi —ipi colors and 
ipi-band magnitudes of stars located more than 18' from 

Assuming they are members of a single stellar population, 
which seems to be the case base d on a lack of intrinsic scatter in 
[Fe/H] of RGB stars (see Section|3.2|). 





14.4 


14.6 


14.8 



Figure 4. Marginal posterior distributions of distance modulus 
at £o = 5° (top) and metallicity Z (bottom). In the bottom panel, 
the solid vertical line shows the mean metallicity Z measured from 
spectroscopy ([Fe/H] = —1.95 di 0.05 dex, [a/Fe] = 0.4 di 0.1 dex), 
while the dashed lines show the uncertainty in the spectroscopic 
estimate of Z. 



Figure 5. Heliocentric distance of the Ophiuchus stream as a 
function of galactic longitude £. The thick solid line shows the 
most probable model, DM(£) = —0.20(£ — 5) + 14.57 mag. To 
illustrate the uncertainty in the most probable model, the thin 
semi-transparent red lines show 200 models sampled from the pos¬ 
terior distribution. The ve rtica l dashed lines show the likely extent 
of the stream (see Section |3.4[ |. The white circles plotted on top 
of the solid line show the positions of 14 confirmed stream mem¬ 
bers, where their distance modulus was calculated using the most 
probable model of DM{£). 


the equatoi p4of th e Ophiuchus stream. 

In Sectioii~ |3.3.2[ we adopted a uniform prior for the 
metallicity content Z in order to test the predictive power 
our dataset. As shown in the bottom panel of Figure]^ 
the peak of the marginal posterior distribution of Z is 
consistent with the mean value of Z estimated from spec¬ 
troscopic data (solid vertical line), and the distribution is 
even narrower than the distribution of Z estimated from 
spectroscopy (dashed vertical lines). This result demon¬ 
strates the predictive power of our data. It shows how 


See Section 3 of Bernard et al.|2014b 


for its definition. 









































Figure 6. The gpi — ipi vs. ipi color-magnitude diagram 
showing the most probable isochrone (yellow thick line) and 200 
isochrones randomly sampled from the stream’s full posterior dis¬ 
tribution (semi-transparent dark red thin lines). The isochrones 
have been shifted to match the distance of the stream at £o = 5°, 
and have been reddened assuming E{B — V) = 0.19 mag (median 
Espd(B — V) at the position of the stream). The grayscale pixels 
show the density distribution of field stars in this diagram (i.e., 
their probability density function). For illustration only, the mag- 
nit udes of observed stars have been corrected for extinction using 
the ISchlegel et al.| l|1998f dust map (to correct for gradients in 
reddening], and tnen again extincted assuming E{B — V) = 0.19 
mag. In addition, the ip^i-band magnitudes of observed stars 
have been corrected for the gradient in distance modulus by adding 
(<? — to) = —0.2(f — 5) mag. Note that the uncertainties in 
color and magnitude also include the uncertainty in extinction. 

a combination of good coverage of the CMD, PSl pho¬ 
tometry, and detailed modeling can provide an accurate 
and precise estimate of the metallicity of single stellar 
populations. 

3.4. Modeling the proper motion and the extent of the 

stream 

The lon gitud e-dependent CMD model we have built 
in Section |3.3[ and the luminosity functions associated 
with the rnodel, allow us to assign a likelihood that a 
star is a member of the Ophiuchus stream, based on the 
star’s galactic longitude €, gpi — ipi color, and ipi-band 
magnitude. The distribution of field stars in the gpi—ipi 
vs. ipi CMD (grayscale pixels in Figure!^, on the other 
hand, enables us to estimate the likelihood that a star is 
associated with the field. As we show in this section, 
these two probability density functions (PDFs), when 
combined with positional and proper motion data, can be 
used to simultaneously trace the extent of the Ophiuchus 
stream and determine its proper motions across the sky. 

In principle, we could measure the proper motion of 
the Ophiuchus stream using the proper motion of its con¬ 
firmed members. However, since our sample of confirmed 
members contains only 14 stars, there is a possibility that 
one or two stars with incorrectly measured proper mo¬ 
tions may bias the results. As an example, stream mem¬ 
ber “rgb4” is clearly an outlier in proper motion as it 
has ~ —24 mas yr“^, while the remaining members 
have mas yr“^. A visual inspection of digitized 

photographic plates has revealed that “rgb4” is blended 


with a neighbor of similar brightness, which affects the 
measured position of the star and its proper motion. 

Fortunately, we do not need to rely only on confirmed 
members and can use a much larger sample of stars in the 
vicinity of the Ophiuchus stream to constrain its proper 
motion and extent. As we deta il be low, we use a prob¬ 
abilistic approach (see Sections 3.1 and 3.31 and model 
the distribution of stars simultaneously in coordinate, 
proper motion, and color-magnitude space as a mixture 
of stream and field (i.e., non-stream) stars. Even though 
we do not a priori know which star is a true member of 
the stream, we assume that as an ensemble, the stream 
stars have certain characteristics which make them dis¬ 
tinguishable from field stars (e.g., common proper mo¬ 
tion, distance, position on the sky and in the CMD), 
and that the scatter in these characteristics is sufficiently 
small to overcome the fact that there are a lot more field 
than stream stars. The narrow width of the stream in 
color-m agnitude (Figure!^ a nd coordinate space (Fig¬ 
ure 1 of Bernard et al.|[20Mb I support this assumption. 
After all, it the stream did not have these characteristics, 
it likely would not have been detected by [Bernard et al.| 
(2014b) in the first place. 

Even though the probabilistic approach we describe 
below uses all of the stars in the vicinity of the Ophi¬ 
uchus stream to constrain the its extent and proper mo¬ 
tion, we expect that most of the signal will come from 
main sequence turn-off (MSTO) stars associated with the 
stream. As shown in Figure the stream’s MSTO is 
bluer than the field population, which means that stars 
in this region of the CMD are much more likely to be 
associated with the stream than with the field popula¬ 
tion. Thus, the statistical weight of such stars will be 
greater than, for example, the weight of stream’s RGB 
stars, which occupy the region of the CMD that is heavily 
dominated by field stars. 

Assuming the stream extends between galactic longi¬ 
tudes imin and imax, the likelihood that a star with 
galactic longitude latitude 6^, proper motions in 
galactic coordinates of Bi,k and color {g — i)k and 
magnitude ik is drawn from the mixture model, is equal 
to 


p(dk|6», 4) = fPstr{dw\dstr,h) + (1 - f)Pfld{dk\6fld, 4), 

( 11 ) 

where dk = {4, bk, Bt,k, Bb,k, [g - *)fc, 4} contains mea¬ 
surements for data point (star) k, and 9 = {0str, Ofid} 
contains parameters that model the distribution of 
stream and field stars, respectively. The parameter / 
specifies the fraction of stars in the stream (out of all 
stars between £min and (-max) and is 0 < / < 1 for all £k 
where imm < £k < £max, otherwise, it is / = 0. 

The likelihood Pstr(dk|6*str4fc) is a product of spatial 
likelihood proper motion likelihood p’fff, and the 
color-magnitude likelihood 

Pstr{dk\9str,£k) = p7tr (&fc I ) 

XpX(W./c,M&.fclX™:4) 

xp<Z{{9-i)kM0Tr7k.bk,EsFD{B-V\£kM))- 

( 12 ) 

The likelihood for field stars, p/zd(dk|0/id4fc)j has the 
same decomposition. 

In galactic coordinates, the distribution of stream stars 













9 


in the latitude direction is modeled with a Gaussian of 
width ab, where the latitude position of the Gaussian 
changes as a quadratic function of the galactic longitude 

=^f{bk\l'{ik),<Jb), (13) 

where iy{ik\A, B,C) = A + B{ik — 4) + C{ik — 4)^ is 
the galactic latitude of the equator of the stream and 

4 = 5°. _ _ 

The spatial distribution of field stars is modeled in a 
similar fashion, except the Gaussian has a width (t(, and 
its latitude position is a linear function of the galactic 
longitude 

pfidih\OfiaJk) = Af{hW'{ik),(Tb), (14) 

where J^'(£fe|y4', S') = A' + B'{£k — io). We use the above 
model for because it is easy to implement, and be¬ 
cause for large ratios of cr(,/crb the Gaussian that models 
the spatial distribution of field stars approximates to a 
plane with respect to the much narrower Gaussian that 
describes the spatial distribution of stream stars. 

At the reference galactic longitude £ = 5°, the stream is 
assumed to have proper motion Jli and Jib, with possible 
gradients in proper motion of ^ and (i.e., gradients 
as a function of galactic longitude). The proper motion 
likelihood of stream stars is then 

J^{pe,k\pe{^k),(Tk) (15) 

X Af {pb^k Ipb {£k ) ; ^/c); 

where /i^(4) = ^(4-4) + Wi and ^^(4) = 

^ (4 — 4) +M6 are the predicted proper motions of the 
stream at galactic longitude ik, and cr(, = fe 

is the quadratic sum of the intrinsic proper motion dis¬ 
persion and the uncertainty in the corresponding proper 
motion of data point k. The purpose of parameter Upm is 
to account for any additional scatter in proper motions 
(e.g., due to unaccounted errors). The proper motion 
likelihood of field stars has the same form (but different 
parameters) as the proper motion likelihood of stream 
stars. 

The likelihood that a star is drawn from the stream’s 
CMD is defined as 

p<Z{{9-i)k,ik\etTr.ik,hk,EsFD{B-V\£k,bk)) = 

cj yW((ff-i)fc|5-b0-(3-i)J 

xA/'( 4 l*,CTiJp (5 - i,i\str)d{g - i)di, 

(16) 


where and are the uncertainty in color and 

magnitude of data point k, and C is a normalization c on- 
stant calculated such that the integral of Equation 
over the considered region of CM space is unity. To ac¬ 
count for the extinction and the gradient in distance, we 
use color {g-iYk = {g-i)k-^A9EsFD{B-V\£k,bk) and 
magnitude*';, = ik-l.682EsFD{B-V\£k,bk)-^^{£k- 
5), where = —0.20 mag deg“^ is the most probable 
gradient in distance modulus (see Table 
In Equation[T^ p{g — i, i\0sYJ) is the PDF of the Ophi- 
uchus stream m the gpi — ipi vs. *pi color-magnitude 


space at galactic longitude 4 = 5°. This PDF was con¬ 
structed by sa mpling isochrones from the stream’s CMD 
model (Section 3.3.3), multiplying them with their lumi¬ 
nosity functions, and then summing them up in a binned 
gpi — ipi vs. ipi CMD. 

The likelihood that a star is drawn from the field CMD 
is calculated as 


pfrAi9-i)k,ik\0t]:g,h,bk) = 

Cj J J^ii9-i)k\9-hcr(g-^)j (17) 

xj\f{ik\i, (r^,Y)p{9 - i, A^lkg, 4, bk)d{g - i)di. 

In Equation [TtI the observed color and magnitude are 
not corrected for extinction or any gradients. 

The PDE of field stars _ (i.e., the CMD), p(g — 

YAb^b'kg^^k,bk) in Equation flTl depends on the galactic 
position. We construct p[g — i,i\9lYYg,£k,bk) by divid¬ 
ing the Ophiuchus region into 1° x 1° spatial pixels that 
overlap by 0.5° in galactic longitude and latitude direc¬ 
tions. Eor each spatial pixel, we bin gpi — ipi colors and 
*Pi"band rnagnitudes of stars located more than 18' from 
the equatoij^ of the stream, and normalize the resulting 
CMD to unity area. An example CMD of field stars is 
shown in Eigure|^as grayscale pixels. 

In total, our model contains 20 parameters. Eor all of 
the parameters, we have adopted uniform priors within 
reasonable bounds. The allowed ranges of model param¬ 
eters were determined by examining positions and proper 
motions of confirmed members and other stars. In addi¬ 
tion to adopted priors, we also require that the parame¬ 
ters satisfy the following constraints: 


1. the spatial width of the stream must be smaller or 
equal than the width of the spatial distribution of 
field stars: Ob < <j'b 

2. the additional scatter in proper motion of stream 
stars must be smaller than the scatter in proper 
motions of field stars: a pm < (^'pm^ 

3. the galactic latitudes of confirmed members {b'jY^^) 
must be within 8ab of the equator of the stream: 
|^(^con/)_^con/| ^ where i*((“"^) is the galac¬ 
tic latitude of the stream at the position of con¬ 
firmed members, and 1 < k < 14. 

As our data set, we use stars brighter than *pi =20 
mag with measured proper motions, and located in a 4 x 4 
deg^ area centered on the Ophiuchus stream. To explore 
the parameter space, we use 200 emcee walkers and ob¬ 
tain convergence after a short burn-in phase of 100 steps. 
The chains are then restarted around the best-fit value 
and evolved for another 2000 steps. The maximum a pos¬ 
terior values, the median and the central 68% confidence 
intervals of model parameters are listed in Table 
We find the stream to be confined between galactic 
longitudes of 3.81° and 5.85° (Figure When com¬ 
bined with the distance of the stream (Figure [^, this 
result implies that the deprojected length of the stream 
is 1.6 ± 0.3 kpc. Thus, the stream is very foreshort¬ 
ened in projection, by a ratio of 6 : 1. The galactic 


See Section 3 of Bernard et al.|2014b 


for its definition 
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Figure 7. The extent of the Ophiuchus stream in galactic coordi¬ 
nates. The gray-scale map shows the probability-weighted number 
density of the Ophiuchus stream, smoothed using a 6^-wide Gaus¬ 
sian filter. The thick solid line shows the most probable model 
(b{i) = 31.37 — 0.80(£ — 5) — 0.15(^ — 5)^ deg) for the equator of the 
stream. To illustrate the uncertainty in the most probable model, 
the thin semi-transparent red lines show 200 models sampled from 
the posterior distribution. The v ertic al dashed lines show the likely 
extent of the stream (see Section [3^. The yellow points show the 
positions of confirmed members, the blue points show candidate 
BHB stars (probability of being stream members > 80%), and the 
arrow indicates the direction of movement of the stream. 



Figure 8. Proper motion of the Ophiuchus stream, inferred from 
the ensemble of likely stream members. The panels show the proper 
motion in galactic longitude {top) and latitude directions {bottom)^ 
as a function of galactic longitude £. The thick solid lines show 
the most probable models {pi(i) = —1.5(^ — 5) — 5.5 mas yr~^, 
fii{£) = 2.0(£ — 5) -|- 2.4 mas yr~^). To illustrate the uncertainty 
in most probable models, the thin semi-transparent red lines show 
200 models sampled from respective posterior distributions. For 
comparison, the semi-transparent gray lines show the proper mo¬ 
tion of field stars. The ver tica l dashed lines show the likely extent 
of the stream (see Section [3.4[ |. 


latitude of the equator of the stream is at bstreami^) = 
31.37 —0.80(£—5) —0.15(£ —5)^ deg, and the width of the 
stream is Cb ~ 6 arcmin (in the galactic latitude direc¬ 
tion) . In direction perpendicular to the stream’s equator, 
the stream is ~ 3.5' wide, which is sim ilar to width of 
^ 3' measured by Bernard et al. (2014b|). 

Using the gpi—ipi color and tp,-band magn itude, and 
given the CMD model of the stream (Section 3.3.31, we 
can evaluate the probability that a star is a member of 
the Ophiuchus stream. We have calculated these prob¬ 
abilities for all of the stars in the vicinity of the Ophi¬ 
uchus stream and have created a probability-weighted 
number density map of the stream, shown as grayscale 
pixels in Figure An inspection of the number den¬ 
sity map did not reveal a significant overdensity of stars 
along the stream that would indicate the presence of a 
progenitor. 

Our data indicate that the proper motion of the stream 


changes as a function of galactic longitude (Figure |^. 
The gradients in proper motion are significant at > ocr 
level, and while their absolute values are similar (~ 2 


mas yr ^ deg ^), the gradients have opposite signs. For 
comparison, the gradients in proper motions of field stars 

M 


are 10 times smaller, 

d)i'h 


~ 0.2 mas yr ^ deg 


T-i 


and = —0.20 mas 


0.1 mas yr ^ deg ^ and 
Overall, the proper mo¬ 
tions of field stars (^^ ^ —6 mas yr“^ and /rj, ~ 0.2 mas 
yr“^) are consistent with apparent motions of a popula¬ 
tion at ^ 8 kpc (due to the motion of the Sun around 
the Galaxy). For comparison, the apparent motion of 
the compact radio source Sgr A* at the Galactic center 
is = -6.38 mas yr"^ 

yr“^ ( Reid fc Brunthaler||2004 ). 

The stream parameters we have obtained so far can 
be used to place a lower limit on the mass of the initial 
population of the Ophiuchus stream. The fraction of 
stars / in the stream between longitudes imin and imax, 
can be converted to the number of stars in the stream, 
Tatars ■ We find that there are Ngtare = 300 ± 30 stars 
brighter than ipi = 20 mag in the Ophiuchus stream. 
If we adopt the luminosity function associated with the 
most probable CMD model of the stream and assume 
Kroupa (19981 initial mass function (not corrected for 
binarity), this number of stars implies that the initial 
population of the Ophiuchus stream had to have a mass 
of at least Minu = (7.0 ± 0.7) x 10^ Mq. 

4. ORBIT OF THE OPHIUCHUS STREAM 

The data and models of the stream obtained in pre¬ 
vious sections now enable us to constrain the orbit of 
the Ophiuchus stream. For this purpose we use galpjj^ 
a package for galacti c dynamics written in Python pro¬ 
gramming language (|Bovy||2015|). 

To make the best use ot the stream constraints and 
their covariances derived so far, we sample the PDF of 
stream constraints with 200 stream samples. Each of 
these samples consists of the line of sight velocity, dis¬ 
tance modulus, galactic position, and proper motion for 
14 stars that uniformly sample the stream in galactic 
longitude from £min to imax^ where these two values are 
drawn for each of the 200 samples from the posterior 
distribution. 

To emulate the width of the stream, we assign an un¬ 
certainty of Ub to positions of data points. To all data 
points we assign a 3% uncertainty in distance, 2 km s“^ 
uncertainty in velocity, and 2 mas yr“^ of uncertainty 
in proper motions. We have verified that our results do 
not change significantly if these uncertainties are changed 
within reason. To convert the observed values into 3D 
positions and velocities, galpy assumes the Sun is lo¬ 
cated 8 kpc from the Galactic center (i?o = 8 kpc), the 
circular velocity at the solar radius is VdrciRo) = 220 km 
s“^, and Sun’s motion in the Galaxy is (-11.1, 244, 7.25) 


km s“^ (Schonrich et al.||2010 Bovy et al.||2012). 

We fit the orbit ot each of the 200 stream samples. 
The orbits are integrated in the default galpy poten - 
tial, called MWPotential2014 (Table 1 of Bovy 2015). 
This potential consists of a bulge modeled as a power- 
law density profile that is exponentially cutoff with a 


http://github.com/j obovy/galpy 
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power-law exponent of -1.8 and a cut-off radius of 1.9 
kpc, a Miyamoto-Nagai disk, and a dark-matter NFW 
halo. MWPotential2014 is consistent with a large vari¬ 
ety of dynamical constraints on the potential of the Milky 
Way, ranging from the bulge to the outer halo. 

The best-fit line of sight velocities, heliocentric dis¬ 
tances, positions, and proper motions predicted by galpy 
orbits for each stream sample are shown as thin semi¬ 
transparent blue lines in Figure The maximum a pos¬ 
terior values, the median and the central 68% confidence 
intervals of orbital parameters are listed in Table 

Overall, the observed and predicted mean values and 
gradients agree within uncertainties. This agreement is 
not trivial. While there is always an orbit that will fit 
a single star in some potential, the same is not true for 
a stream of stars. For example, given the observed gra¬ 
dients and mean values in proper motion, distance, and 
position, the observed gradient in line of sight velocity 
has to be positive, otherwise, there is a strong discrep¬ 
ancy with the velocity predicted by the most probable 
orbit. Similarly, the observed gradient in distance mod¬ 
ulus has to have a negative sign, otherwise a plausible 
orbit fit cannot be achieved. 

The most noticeable disagreement is between observed 
and predicted proper motions (bottom right panel of Fig¬ 
ure^, with the observed proper motion in the longitude 
direction at £o = 5° iWe) having a —2.2 mas yr“^ off¬ 
set with respect to the proper motion predicted by the 
galpy orbit ht. The result of such observed proper mo¬ 
tion (i.e., Jle = —5.5 mas yr“^) is that the velocity vectors 
of stream stars do not align with the extent of the stream 
in the galactocentric X vs. Y (and Y vs. Z) plane (see 
Figure 101. The expected behavior would be for the ve¬ 


locity vectors to be aligned with the extent of the stream 
(i.e., the stream gets longer in the direction it is moving), 
which would happen for Jli = —7.7 mas yr~^ . 

As we have already stated in Section |2.4[ candidate 
QSOs do not show any statistically signihcant proper 
motion = 0.3 ± 0.2 mas yr“^), and galaxies do 

not show any bulk motion either. Thus, we have no 
indication that faulty proper motions are the cause of this 
inconsistency between the velocity vector of the stream 
and its extent. 

We have repeated orbit htting after adding a —2.2 mas 
yr“^ offset to proper motions in the longitude direction, 
and have found that our results do not change. This was 
expected since we assume a 2 mas yr“^ uncertainty in 
proper moti ons when fitting orbits with galpy. 

Figure [^illustrates the orbit of the Ophiuchus stream 
in the past 350 Myr. We hnd that the stream has 
a relatively short orbital period of 346^7^ Myr, and a 
fairly eccentric orbit (e = O.Gb^Q giJ), with a pericenter of 

3.57lg;o6 kpc and an apocenter of 16.8 ^q 4 kpc. About 
10 Myr ago, the stream passed through its pericenter 
and now it is moving away from the Galactic plane and 
towards the Galactic center. 

The above uncertainties in orbital parameters only ac¬ 
count for the uncertainties in position and velocity of the 
stream, and do not account for the uncertainty in the dis¬ 
tance of the Sun from the Galactic center (i?o), and the 
circular velocity at the solar radius (vdrciRo))■ To deter¬ 
mine how the orbital parameters change as a function of 


i?o and VcirciRo), we fit the orbit of the stream assuming 
a Milky-Way-like potent ial fit to dyn amical data as de¬ 
scribed in Section 3.5 of Bovy (2015), but assuming val¬ 
ues of (8.5 kpc, 220 km s“^), (8.0 kpc, 235 km s“^), and 
(8.5 kpc, 235 km s“^) for {RQ,Vcirc{Ro))- When fitting 
these orbits, we model the stream using the maximum a 
posterior values listed in Table [I] 

We find that the best ortat fit is obtained for 
(Ro^VcirciRo)) = (8.0,220) (i.e., the default galpy val¬ 
ues). With respect to fiducial periods (i.e., those ob¬ 
tained assuming Rq = 8.0 kpc and VdrciRo) = 220 km 
s“^), modifying Rq and Vcirc{Ro) changes periods of the 
Ophiuchus stream between -|-2 and —75 Myr. Other or¬ 
bital parameters do not change appreciably. 

5. TIME OF DISRUPTION 

In the Introduction, we said that the short length of 
the Ophiuchus stream suggests that its progenitor must 
have been disr upted fairly recently. As we have shown 
in Section |3^ part of the reason the stream is so short 
in projection, is the viewing angle-we are observing the 
stream almost end-on. 

Even when the projection effects are taken into ac¬ 
count, the deprojected length is still fairly short, only 1.6 
kpc. For compariso n, the second short est stellar stream 
is the Pisces stream (iMartin et al.|2013l 

3 , 11201 ^ 


also known as the 
Triangulum stream, [Bonaca et al.|2U12 ) with a length of 
5.5 kpc. Therefore, the length of the stream still sug¬ 
gests that the stream formed recently, that is, it suggests 
that the progenitor was recently disrupted. 

As the progenitor of the Ophiuchus stream orbited the 
Galaxy, it would have experienced the tidal force of the 
Galactic potential. This force can strip stars from the 
progenitor and it could have been strong enough to com¬ 
pletely disrupt the progenitor. 

In order to examine the influence of the tidal force, we 
have calculated its magnitude as a function of time for 
the most probable orbit of the Ophiuchus stream. The 
magnitude of the tidal force was calculated by finding 
the largest eigenvalue of the following matrix 


d^<I> d- 




\ dZdR dZ^ J 


(18) 


where 4) is value of the galpy MWPotential2014 poten¬ 
tial at the position of the progenitor, and R and Z are 
coordinates in the cylindr ical galactocentric system. The 
result is shown in Figure 

We find that the tidalTorce is the strongest during 
pericenter-l-disk passages, and that the progenitor of the 
stream could have been disrupted during one of those 
passages. To find if the progenitor could plausibly have 
been disrupted during one of these passages, we can use 

galpy- 

Given an orbit, the time of disruption tdis, and the 
velocity dispersion of the progenitor cr„, galpy can gen¬ 
erate a moc k stream using the modeling framework of 
Bovy (2014). For a fixed (jy, the disruption time, tdis, 
IS proportional to the s trea m’s length (i.e., older streams 
are longer). In Section 3.1 we measured s = 0.4 km s“^ 
as the median velocity dispersion of the stream. We find 
that by setting cr„ to the same value, the galpy mock 
stream has a velocity dispersion s « 0.4 km s“^ = s. We 
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Figure 9. This plot compares line of sight velocities (top left), distances {top right), positions {bottom left), and proper motions {bottom 
right) calculated by galpy (thin blue lines) with models derived from observations (thin red lines). The observed and calculated values 
are consistent within uncertainties, with the exception of the proper motions, where the model cannot match the apparent gradient of the 
proper motions along the stream. 



Figure 10. This plot illustrates the misalignment between the 
observed velocity vectors of stream stars {arrows) and the extent 
of the stream {solid line) in the galactocentric Cartesian X vs. Y 
plane (i.e., a top-down view of the Galactic plane). In this coordi¬ 
nate system, the Sun is at {X, Y, Z) = (8, 0, 0) kpc and the y-axis 
is positive toward galactic longitude I = 270°. 


further find that tdis ^ 170 Myr provides a good match 
between the observed and mock streams (see Figure 13). 
While we do not perform an exhaustive search of the 
(o'vjtdis) parameter space, very different values of or 
tdis produce values of s or a total length that are not 
consistent with observations. This result, and the fact 
that there was a disk+pericenter passage 240 Myr ago, 
strongly suggest that the stream formed about 240 Myr 
ago (i.e., that th e progenito r was d isrupted at that time). 

As shown by Johnston (1998), the velocity disper¬ 
sion of the progenitor scales with the mass of the pro- 
genitor as cr^ oc In galpy, the model used for 

stream generation was calibrated using a progenitor of 
mass Mdyn = 2 x 10^ Mq. The stream formed by the 
tidal disruption of t his progenit or could be modeled with 
cr„ = 0.365 km s“^ (Bovy|2014). Using the above scaling 



X [kpc] 

Figure 11. The orbit of the Ophiuchus stream in the past 350 
Myr (about one orbital period), shown in a right-handed galac¬ 
tocentric Cartesian coordinate system. In this coordinate system, 
the Sun is at {X, Y, Z) = (8, 0, 0) kpc and the y-axis is positive 
toward galactic longitude I = 270°. Note the pericenter passage at 
t ^ —240 Myr (solid square). Near this point in time, the stream 
was also passing through the disk (Z ~ 0 kpc, see the top left panel) 
and w as e xperiencing strong tidal forces due to disk shocking (see 
Figure [T^. 

relation and = 0.4 km s“^, we find that the progenitor 
of the Ophiuchus stream had a mass of ^ 2 x 10^ 

It is important to note that galpy creates only the 
stream, and that stars associated with the progenitor 
are not part of the mock stream (i.e., are not shown in 
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Figure 12. The tidal force acting on the Ophiuchus stream, 
normalized to the maximum tidal force in the past 2 Gyr. The 
narrow p eaks correspond to p assages through the disk (i.e., disk 
shocking, |Ostriker et al.||1972|| and the broader peaks correspond 
to passages tnrougn tne pericenter. 



Figure 13. A comparison of the observed number density map 
of the stream (top), a map created from a mock stream generated 
using galpy (middle), and a map created from a Af-body stream 
generated using NEMO. In all panels, the solid line shows the most 
probable position of the stream and the dashed lines illustrate its 
1(7 width. The mock galpy stream was generated assuming time 
of disruption = 170 Myr, and velocity dispersion = 0.4 km 
s“^. Note a good agreement between the length and width of the 
observed, galpy and NEMO streams. 


the middle panel of Figure 13). This is the reason an 
overdensity of stars (i.e., the progenitor) is not visible in 
the mock galpy stream. 

To create a more realistic stream that includes a pro- 
genitor, we use the gyrf alcON code (|D ehnen|2000l [2002 1 
in the NEMO toolkit (|Teub en||1995|). We set up the 
progenitor as a King cluster (King |l966 ) with a mass of 
1X 10^ Mq, tidal radius of 94 pc and a ratio of the central 
potential to the velocity dispersion squared of 2.0. The 
cluster is sampled using 20,000 particles and is evolved 
for ~ 365 Myr in the MWPotential2014 potential. The 
initial conditions are those at the t ~ —365 Myr apocen- 
ter of the most probable orbit of the Ophiuchus stream. 

The number density map of the A^-body stream created 
using gyrfalcON is shown in the bottom panel of Fig¬ 
ure The width and the length of the IV-body stream 
match the observed stream fairly well. For comparison, 
if the cluster is evolved starting from the apocenter at 
t ^ —870 Myr, the resulting stream is much longer and 
inconsistent with observations. Based on this more real¬ 
istic simulation, we conclude that the Ophiuchus stream 
likely formed about 240 Myr ago and that its progeni¬ 
tor was fully disrupted during a single disk-bpericenter 
passage. 

6 . CONCLUSIONS AND SUMMARY 

In this paper, we have presented follow-up spec¬ 
troscopy and an astrometric and photometric analysis of 
the Ophiuchu s stellar stream i n the M ilky Way, recently 
discovered by Bernard et al. (2014b) in PSl data. We 
have been able to put together a comprehensive, empir¬ 
ical description of the Ophiuchus stream in phase space: 
we succeeded in determining the mean phase-space co¬ 
ordinates in all six dimensions, along with the gradients 
of those coordinates along the stream (see Table for a 
summary of stream parameters). 

Overall, phase-space data along the stream can be well 
matched by an orbit in a fiducial Milky Way potential: 
Ophiuchus is truly a thin and long (~ 1.6 kpc) stellar 
stream, 50 times longer than wide, that appears 6 : 1 
foreshortened in projection; it is on a highly inclined or¬ 
bit with only a 350 Myr orbital period; it is receding 
from us at nearly 300 km s“^ and has just passed its 
pericenter at ~ 3 kpc from the Galactic center. This 
makes Ophiuchus the innermost stellar stream known in 
our Galaxy. It is also the only known kinematically cold 
stellar stream to be seen nearly end-on. 

The homogeneously metal-poor ([Fe/H] = —2.0 dex), 
a-enhanced ([a/Fe] ^ 0.4 dex) and old stellar population 
(^ 12 Gyr old), and the small line of sight velocity disper¬ 
sion we found (< 1 km s“^), strongly confirm the notion 
that the progenitor of the stream must have been a glob¬ 
ular cluster. If the detected part of the stream encom¬ 
passes most of the progenitor’s stars, then the progeni¬ 
tor’s tidal radius was ^ 90 pc and its mass was ~ 2 x 10"^ 
Mq (certainly greater than ~ 7 x 10^ Mir,) . In this re¬ 
spect, the Oph i uchus and the GD-1 s tream (jGrillmair &| 
Dionatos] 2006 [Koposov et al.j 2010) can be considered 
identical twins, as they have the same metallicity, the 
same mass, and no detectable progenitors. 

Our analysis, however, leaves a number of questions 
open. First, the most probable orbit in the hducial po¬ 
tential is not quite able to match the proper motions and 
their gradients along the stream. A thorough exploration 
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whether there are axiymmetric or non-axisymmetric po¬ 
tentials that might be able to remedy this tension re¬ 
mains to be done. Alternatively, this discrepancy may 
indicate a problem with the proper motion data. While 
we have done our best to obtain good proper motions, we 
cannot fully dismiss this possibility. However, we expect 
that the prope r motions provided by the the ongoing 
GAIA mission (Perryman et al. 20011 will resolve this 
discrepancy in the near future. 

Second, the present analysis does not yet use the 
stream phase-space data to provide new constraints on 
the Galactic potential. In principle, the Ophiuchus 
stream can provide constraints on the Galactic poten¬ 
tial at about 4 kpc above the Galactic center, a location 
where few other constraints exist. As the top left panel 
of Figure Pshows, the line of sight velocity of the stream 
is predicts to decrease to about 282 km s“^ at galactic 
longitude Z = 4°. By measuring line of sight velocities 
of stars at this position, we can identify members of the 
Ophiuchus stream and test the assumed potential. Sim¬ 
ilarly, the potential may also be tested by identifying 
stream members at I > 5.8°, where the stream is pre¬ 
dicted to curve toward the galactic latitude b = 30.5° 
(see the blue lines in the bottom left panel of Figure]^. 

And finally, our analysis suggests that the progeni¬ 
tor of the stream was fully disrupted during a single 
disk-l-pericenter passage about 240 Myr ago. As IV-body 
simulations show, this disruption was strong enough to 
smooth out the distribution of stars along the stream 
and effectively erase all evidence of the progenitor. If 
this scenario is correct, then the answer to the question 
“Where is the progenitor of the Ophiuchus stream?” is 
fairly simple-the Ophiuchus stream is all that is left of 
the progenitor. 

While the above scenario seems to answer one ques¬ 
tion, it creates another one. The fact that the ?V-body 
stream ends up being too long if the cluster is evolved 
for more than ~ 400 Myr, suggests that the progenitor 
could not have been on the current orbit for more than 
^ 400 Myr. If that is true, how did the progenitor end 
up on the current orbit and what was its original orbit? 
We expect that detailed A^-body simulations that include 
the interactions with the Galactic bar will provide a more 
definitive answer to these questions and plan to pursue 
this approach. 

Our finding that the progenitor could not have been 
on the current orbit for more than ^ 400 Myr is based 
on the comparison of the length of the observed and the 
A^-body stream. However, what if the stream is actually 
much longer, but is simply not observed as such in cur¬ 
rent data? This could happen if, for example, the stream 
suddenly fans out and thus its surface brightness drops 
below our detection limit. The stream may fan out due 


panel of Figure 3 by 

Bonaca et al. 

2014), due to being 

on a chaotic orbit (bottom panels o 

Figure 11 by Fardal 

et al. 

2015 also [Price-Whelan et al. 

2015), or due to 

being 

in a triaxiai potential (l:^'igure 4 by Fearson et al. 


20151. A longer stream would imply that the progenitor 


has been undergoing disruption for a longer time, which 
would make the observed orbit more plausible as the pro¬ 
genitor’s original orbit (i.e., a change in orbit would not 
be necessary). 

The true extent of the stream may be constrained by 


identifying stream members along the predicted extent of 
the stream or in a fan-out pattern, via line of sight veloci¬ 
ties. As Figure[2 shows, the stream’s high velocity makes 
the separation ra member stars from field stars an easy 
task. We have already started a follow-up spectroscopic 
program with the goal of identifying additional stream 
members, and hope to better constrain the length, orbit, 
and possible fanning-out of the Ophiuchus stream in the 
near future. 
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Table 4 

Ophiuchus Stream Member Stars 


Name 

R.A. 

(deg) 

Decl. 

(deg) 

gpi 

(mag) 

rpi 

(mag) 

ipi 

(mag) 

Zpi 

(mag) 

ypi 

(mag) 

'^los 

(km s“^) 

DM 

(mag) 

(mas yr“^) 

fJ'b 

(mas yr“^) 

bhbl 

241.52271 

-7.01555 

16.05 ±0.02 

16.11 ± 0.02 

16.22 ±0.01 

16.25 ± 0.02 

16.25 ±0.02 

286.7 ± 1.8 


-4.1 ± 2.1 

2.6 ±2.1 

bhb2 

241.49994 

-7.03409 

16.02 ±0.02 

16.09 ±0.02 

16.21 ± 0.02 

16.25 ± 0.02 

16.23 ±0.02 

285.3 ± 1.9 

14 73+0 '^® 

<0_Q Qg 

-4.4 ± 2.2 

2.1 ± 2.2 

bhb3 

242.13551 

-6.87785 

15.96 ±0.02 

15.99 ±0.01 

16.11 ± 0.02 

16.16 ±0.02 

16.13 ±0.02 

290.0 ± 1.8 

14 61"^®'®^ 

Q5 

-5.2 ±2.0 

2.4 ± 2.0 

bhb4 

241.94714 

-6.88995 

16.22 ±0.02 

16.32 ±0.01 

16.48 ± 0.02 

16.54 ±0.02 

16.56 ±0.02 

291.3 ±2.2 


-5.1 ± 2.2 

3.3 ±2.1 

bhb6 

242.33018 

-6.84405 

15.67 ± 0.01 

15.64 ± 0.01 

15.59 ±0.01 

15.60 ±0.02 

15.54 ± 0.02 

290.8 ± 1.5 

i4.00_Q Q5 

-5.6 ± 1.3 

5.2 ± 1.4 

bhb7 

242.91469 

-6.69329 

15.66 ±0.01 

15.56 ±0.02 

15.57 ± 0.01 

15.54 ±0.02 

15.50 ± 0.02 

289.8 ± 1.5 

14 47+0 05 

-7.4 ± 1.5 

7.1 ± 1.7 

bhb6 








289.8 ±0.8 




rgbl 

241.51689 

-6.98511 

17.71 ± 0.02 

17.18 ±0.02 

16.91 ± 0.02 

16.78 ±0.02 

16.71 ± 0.02 

286.0 ±0.8 

-L++^_0.06 

-6.2 ± 1.5 

0.7 ± 1.4 

rgb2 

241.94649 

-6.86113 

17.34 ± 0.02 

16.74 ± 0.02 

16.46 ± 0.02 

16.32 ± 0.02 

16.25 ±0.02 

286.7 ± 0.6 

14 fi4+0 05 

i4.D4_Q Q5 

-4.5 ± 1.4 

1.1 ± 1.4 

rgbS 

241.96089 

-6.89873 

17.64 ± 0.02 

17.08 ±0.02 

16.83 ±0.02 

16.68 ±0.02 

16.62 ±0.02 

287.5 ±0.7 

14 fi4+0 05 

i4.D4_Q Qg 

-6.5 ± 1.5 

3.0 ± 1.4 

rgb4 

242.26139 

-6.90190 

17.05 ±0.02 

16.45 ±0.02 

16.17 ± 0.01 

16.02 ± 0.02 

15.93 ±0.02 

288.8 ±0.5 

14 60 

i4.DU_Q Q5 

-23.7 ± 1.4 

-14.8 ± 1.5 

rgb5 

242.14832 

-6.79765 

17.79 ±0.02 

17.23 ±0.02 

16.96 ±0.02 

16.82 ± 0.02 

16.74 ± 0.02 

288.0 ±0.9 

i4.DU_Q Q5 

-6.0 ± 1.3 

2.4 ± 1.4 

sgbl 

241.95962 

-6.86881 

18.90 ± 0.02 

18.53 ±0.02 

18.38 ±0.02 

18.31 ± 0.02 

18.28 ±0.02 

289.4 ± 2.2 

14 «3+0 05 
i4.D<5_Q 05 

-5.8 ±2.0 

1.2 ±2.1 

mstol 

242.02040 

-6.84122 

19.22 ±0.02 

18.89 ±0.02 

18.74 ±0.02 

18.69 ±0.02 

18.62 ±0.03 

291.8 ±2.2 

i4.DZ_0.05 

-6.9 ±2.4 

4.4 ± 2.4 

msto2 

242.18360 

-6.84056 

19.10 ± 0.02 

18.70 ±0.02 

18.54 ±0.02 

18.45 ± 0.02 

18.42 ±0.03 

286.4 ± 2.6 

i4.DU_o.o5 

-3.4 ± 2.3 

2.3 ±2.2 


Note. — The horizontal line separates stars observed by DEIMOS and Hectochelle. The name indicates the likely evolutionary stage inferred from isochrone fitting. 
The PSl photometry is not corrected for extinction. The uncertainty in vios includes the uncertainty from cross-correlation/fitting and the uncertainty in the zero-point of 
wavelength calibration. The DM indicates the average DM of the stream at the position of the star, and the uncertainties are 68% confidence intervals. The proper motion 
in the galactic longitude direction, /x^, includes the cos 6 term. 



